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ABSTRACT 

RX J1856.5— 3754, a bright soft X-ray source believed to be the nearest thermally emitting neutron 
star, has commanded and continues to command intense interest from X-ray missions. One of the 
main goals is to determine the radius of this neutron star. An integral part of the determination is an 
accurate parallax. Walter (2001) analyzed Hubble Space Telescope (HST) data and derived a parallax, 
7r = 16.5 ± 2.3 mas. Combining this distance with the angular radius derived from blackbody fits to 
observations of RX J1856.5-3754 with ROSAT, EUVE, HST, Pons et al. (2002) derived an observed 
radius ("radiation radius"), R^ = 7km. This value is smaller than the radii calculated from all proposed 
equations-of-state (EOS) of dense baryonic matter (Haensel 2001). Here, we have analyzed the same 
HST data and find ir — 7± 2 mas. We have verified our result using a number of different, independent 
techniques, and find the result to be robust. The implied radius of RX J1856.5— 3754 is Roo = 15 ±6 km, 
falling squarely in the range of radii, 12-16 km, expected from calculations of neutron star structure 
for different equations of state. The new distance also implies a smaller age for RX J1856.5— 3754 of 
0.4 Myr, based on its association with the Upper Sco OB association. 

Subject headings: astrometry — pulsars: individual (RX J1856.5— 3754) — stars: neutron 



1. INTRODUCTION 

The ROSAT all-sky survey identified six neutron stars 
that are radio-quiet but bright in the soft X-ray band. 
These sources, unlike the well studied radio pulsars, lack 
significant non-thermal emission and are thus excellent 
candidates for X-ray spectroscopic studies of the atmo- 
spheres of neutron stars (for reviews, see Motch 2000; 
Treves et al. 2000). 

The brightest of these sources is RX J1856.5— 3754 
(Walter, Wolk, & Neuhauser 1996). A faint, blue opti- 
cal counterpart was identified from Hubble Space Telescope 
(HST) data (Walter & Matthews 1997). RX J1856. 5-3754 
has been intensively studied by most major facilities, espe- 
cially ROSAT, EUVE, ASCA and HST. The broad-band 
data can be well fitted by thermal emission from a neu- 
tron star, which has resulted in the determination of the 
effective temperature, kT c g w 50 eV, and angular radius 
of the neutron star, Roo/d ~ 0.11 km pc _1 (Pavlov et al. 
1996; Pons et al. 2002). 

Walter (2001, hereafter W01) used HST data span- 
ning three years to measure the astrometric parameters 
of RX J1856.5— 3754, finding a parallax of 16.5 ± 2.3 mas 
and a proper motion of 332 ± 1 mas yr _1 at a posi- 
tion angle of 100.3° ± 0.1° (/i Q = 326.7 ± 0.8 mas yr _1 , 
= —59.1 ± 0.7 mas yr _1 ). 

Combining this parallax with the broad-band modeling 
yields a radiation radius of Roo ~ 7 km. For the canonical 
mass of a neutron star, 1.4 Mq this radius is smaller than 
the minimum radiation radius of i? 1 " 111 = 10.7 km allowed 



by General Relativity (Lattimer & Prakash 2000). For the 
radius to exceed i?™ ln the mass has to be less than 1 Mq 
(Pons et al. 2002). 

The importance of RX J1856.5— 3754 as a laboratory for 
dense matter physics has motivated deep observations by 
Chandra (Burwitz et al. 2001), XMM and the Very Large 
Telescope (VLT). Results from the first 50-ksec Chandra 
observations 1 can be found in Burwitz et al. (2001); the 
blackbody fits are similar to those of Pons et al. (2002). 
Using the VLT, van Kerkwijk & Kulkarni (2001b) dis- 
covered an unusual Ha nebula around RX J1856.5— 3754, 
from which they infer properties of RX J1856.5— 3754's 
energetics and emission characteristics (van Kerkwijk & 
Kulkarni 2001a). 

Parallax or distance is essential to obtaining the radius, 
the key physical parameter (since it now appears that X- 
ray and optical observations yield reliable values for T e g 
and the angular radius). In view of the perplexing radius 
inferred from the parallax measurement of W01 we under- 
took an analysis of the publicly available HST data. Here 
we present a detailed description of our analysis followed 
by our measurement of the parallax. 

2. OBSERVATIONS, ANALYSIS & RESULTS 

We analyzed the publicly available fiS"77WFPC2 obser- 
vations described by W01 and Walter & Matthews (1997); 
see Table 1 for a log of the observations. As noted by W01, 

1 In early October, 2001, Chandra observed RX J1856.5-3754 for 
an additional 450-ksec under the aegis of the Director's discretionary 
program. 
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the observations took place near the times of parallactic 
maximum (3 October and 30 March). RX J1856. 5-3754 
is always on the Planetary Camera (PC) detector so we 
only analyzed those data. 

2.1. Relative Astrometry 

We used the effective point-spread function (ePSF) fit- 
ting technique to perform the astrometry, as described by 
Anderson & King (2000, hereafter AK00). We did not have 
sufficient numbers of stars to derive our own ePSF for each 
data set, so we used a previously determined ePSF (from 
archival data) for the #ST/WFPC2 F555W filter. While 
the parallax data were taken with the F606W filter, we feel 
that using the F555W ePSF was appropriate, as it was of 
superior quality to the F606W ePSF that we have (also de- 
rived from other data). The difference in ePSFs should not 
bias the data, as the wavelength dependence of the ePSF 
is not strong, especially across the w 50 nm difference in 
effective wavelength between the filters. Furthermore, the 
blue color of RX J1856. 5-3754 brings its ePSF closer to 
the F555W ePSFs of normal stars. In any case, we also 
performed the analysis with the F606W ePSF. Since the 
latter ePSFs were of inferior quality (owing to a less ideal 
data set), we obtained larger errors, but the results were 
entirely consistent with those obtained using the F555W 
ePSF. 

We fit the F555W ePSF to the raw images, uncorrected 
for dithering or cosmic rays. For each epoch, we used a x 2 - 
minimization, as described by AK00, to derive a position 
for each star in each of the raw images. We corrected 
this position for the 34th-row anomaly 2 (Anderson & King 
1999) and geometric distortions in the PC detector using 
new coefficients (Anderson 2002, in preparation). For each 
epoch, this yielded four (or eight) positions for each star. 

We then solved for the shifts (due to dithering) between 
the four (or eight) images in a given epoch; these shifts 
are given in Table 2. We rejected sources which had sig- 
nificantly higher residuals than other sources of their mag- 
nitude (see Figure 1), ascribable to the source being ex- 
tended or saturated. For the remaining sources we used an 
iterative cr-clipping (with threshold at 2.5cr) to reject out- 
lier position determinations. Following this the remaining 
position measurements were used to derive the average po- 
sition for each source. The final distortion-corrected source 
positions in image (x, y) coordinates as well as the number 
of accepted measurements are listed in Table 3. 

In the limit of a large number of independent observa- 
tions that are well dithered, the measured rms would have 
yielded reliable uncertainties. However, we are limited 
by the meager number of available frames at each epoch. 
With a maximum of four (or even eight) measurements of 
each position we were not able to search for systematic er- 
rors. Furthermore, the data were taken with non-optimal 
dithering 3 , with the fractional pixel portions of the dither 
often repeating from one image to the next (Table 2). This 
increases the likelihood of systematic pixel-phase errors. 

2 It is a common error to apply the 34th-row correction after shifting 
and combining the data. This procedure is incorrect, as the 34th-row 
correction should be applied to the raw image coordinates and not 
those that have been shifted and rebinned. 

3 For astrometric purposes, a 2 X 2 dithering grid is minimal for 
optimal removal of pixel-phase errors; a 15-point dither pattern is 
even better. See AK00. 



We therefore adopted a semi-empirical approach for the 
measurement uncertainties along the following lines. 
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FlG. 1. — 1-D position uncertainty (pixels) vs. instrumental mag- 
nitude, defined by m = — 2.5 log 10 (DN) within a 5 X 5 pixel region 
in a single exposure. The x's are for the raw uncertainties in the 
x positions of the stars used to register the different epochs, the 
triangles for the raw y uncertainties, and the solid line is a relation 
determined from ~ 5000 well-observed stars in other data sets. The 
sources here are those that are not saturated and were used in the 
analysis; see Tables 4 and 3. The data from this paper (x's and 
triangles) generally follow the trend denned by the line (also given 
in Eqn. 1), but there is considerable spread due to the small number 
of measurements (3 or 4) used to construct each uncertainty. 

As the first approximation of the uncertainty for each 
position, we take the rms variation between the positions 
used to construct the average. As the next level, we used 
the expected precision in the positions as a function of 
the signal-to-noise ratio (SNR) of stars. To this end, we 
utilized an astrometric database that one of us (J. A.) has 
built up over the last several years. In particular, we used 
18 well-dithered PC images that were obtained under sim- 
ilar conditions (filter, background, crowding) to those dis- 
cussed here. In Figure 1 we display a fit (obtained from 
the measurements of about 5000 stars) to the astrometric 
precision as a function of the SNR of stars. The relation 
from Figure 1 is reasonably well fit by 
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where m = — 2.5 log 10 (DN) (within a 5 x 5 pixel area) is 
the magnitude measured in a single exposure and a is the 
1-D position uncertainty. For N well-dithcrcd exposures, 
the uncertainty is a/V~N, as expected (AK00). 

We see that the raw uncertainties generally follow the 
expected trend, but that there is substantial scatter. This 
is not surprising, given that we may be computing the 
uncertainties from four or fewer than four measurements. 
Therefore, in our analysis we use the maximum of the em- 
pirically determined uncertainty for an individual star and 
the uncertainty from the relation in Figure 1 corrected to 
the appropriate magnitude; the stars that have uncertain- 
ties larger than those inferred from this relation do so due 
to cosmic rays or proximity to bright sources. 

The above analysis gives us reliable and accurate mea- 
surements of stellar positions, but while these stars are in 
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Table 1 



WFPC2 Observation Summary 



Epoch 


^exp 


«dith cr a 


Sky Level 


Nominal 


APA b 


Scale 




(UT) 


(s) 




(DN s- 1 ) 


PA (deg) 


(deg) 


(mas pixel -1 ) 


pixels 


1996 Oct 6 


4800 


4 


0.006 


129.38 


0.012(2) 


45.5936(8) 


408.004(5) 428.143(5) 


1999 Mar 30 


7200 


8 


0.006 


-51.75 


0.0212(7) 


45.6028(5) 


429.221(3) 437.824(3) 


1999 Sep 16 d 


5191 


4 


0.005 


124.97 





45.5938 


417.948 436.803 



a The dither pattern was in sets of four images, with (Ax,Ay) =(0,0), (0,0), (5.5,5.5), and (5.5,5.5) pixels. 

b Defined as the difference between the fitted position angle and the nominal position angle (from the image header 
and Table 1). 

c x and y pixel coordinates of the pointing center, which has a = -18 h 56 m 35?374 and <5 = -37°54'31'.'71; sec §2.3. 
d The values for this epoch were assumed to be correct. 



Note. — See Eqn. 2 for the sense of the transformation. 



Table 2 



Fitted Dither Positions 



Epoch 


Image 


Ax 


Ay 




Number 


(pixels) 


1996 Oct 6 


1 


0.000 


0.000 




2 


-0.049 


-0.112 




3 


5.451 


5.417 




4 


5.416 


5.270 


1999 Mar 30 


1 


0.000 


0.000 




2 


0.022 


0.115 




3 


5.456 


5.455 




4 


5.465 


5.476 




5 


-0.066 


-0.142 




6 


-0.097 


-0.064 




7 


5.368 


5.628 




8 


5.399 


5.532 


1999 Sep 16 


1 


0.000 


0.000 




2 


-0.159 


-0.027 




3 


5.515 


5.542 




4 


5.386 


5.422 
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the background relative to RX J1856.5— 3754 they can still 
have their own motions that will bias our determinations. 
Therefore, to have some idea of the absolute motion of 
the stars in the image, we included in the data measure- 
ments of the positions of two slightly extended sources 
(presumed to be galaxies) present on the HST images (see 
Table 4). As these sources are non-stellar, we could not 
use the AKOO technique to measure their positions. In- 
stead we fitted Gaussian profiles and then applied the same 
distortion corrections as with the other technique. Gaus- 
sian fitting is inherently less accurate than ePSF fitting 
(AKOO), but the errors are important primarily for sources 
that are undersampled by the WFPC2 pixels (i.e. where 
pixel-phase errors are important). The galaxies were rea- 
sonably well resolved (FWHM w 3.2 pixels for source #20, 
FWHM ss 3.6 pixels for source #104), so they should not 
suffer from systematic errors related to undersampling. 

2.2. Registration of Epochs 

To determine the transformation of the background 
sources (all sources except for RX J1856.5— 3754) between 
epochs, we proceeded iteratively. Our basic input data set 
was the 27 stellar sources that had consistent measure- 
ments in each epoch plus the two galaxies (§2.1), given 
in Table 3. The faintest of these sources were as faint as 
RX J1856.5— 3754, and the brightest were w 160 times as 
bright as RX J1856.5— 3754 (the brightest non-saturated 
sources on the WFPC2 images). First, we set the fidu- 
cial positions of the sources to their measured positions at 
epoch 1999.7. We chose this epoch as the effects of par- 
allax between it and epoch 1996.8 are minimized (due to 
similar parallactic angles) and the effects of proper motion 
between it and epoch 1999.3 are minimized (due to close 
proximity in time), thus yielding the best matches to the 
other epochs given the limited information that we have. 
We assumed that the position angle, scale, and pointing 
center of this fiducial epoch are known. The pointing cen- 
ter has no impact on the results, and is simply a dummy 
parameter. The position angle and scale are known to 
reasonable precision (< 0.1° for the position angle, and 
< 0.1% for the scale). For our nominal values, we take the 
values from the HST image headers (listed in Tables 1). 
If these values are wrong at certain levels, they would in- 
troduce errors on those levels into our measurements (i.e. 
if the fiducial scale were wrong by 0.1%, our proper mo- 
tions would be wrong by the same amount) but these are 
systematic effects that are well below the measurement 
precision of the current data. 

For our transformation between the epochs, we assumed 
a simple rotation, scale, and offset. This is valid if the dis- 
tortion correction (§2.1) removed all nonlinear terms. The 
transformation equation between the measured positions 
of star i at epoch j, {xij,yij), to the celestial position 
{Aa it j, ASij) is 

Aajj = — scale, ({xij — Xoj) cosPAj — [y^j — y j) sin PA.,) 

ASij = scale-,- ((x itj - x ,j) sinPA., + (y lJ - y 0tj ) cosPAj)(2) 

where scale, is the plate scale (arcsec per pixel), (xo,j, yoj) 
are the offsets, and PAj is the position angle of epoch j, 
all of which are given in Table 1. 

We performed a x 2 fit between the fiducial positions and 
the positions at the three measured epochs, varying the 



scale, position angle, and offsets of the other epochs. This 
fit gave relatively large % 2 values, due to proper motion 
between the epochs. 

We then fit for updated fiducial positions and proper 
motions of the background sources based on a linear least- 
squares technique (for the galaxies, the proper motion was 
forced to be 0). These positions and proper motions were 
used to re-calculatc the expected positions in the non- 
reference epochs, which dramatically lowered the % 2 val- 
ues. We iterated this procedure (fitting for the transfor- 
mation between the epochs, then fitting for the individual 
positions and proper motions) making sure that the solu- 
tion was converging. After five iterations, \ 2 changed by 
0.2; we considered the solution to have converged. The 
results of the fitting (the fiducial positions and proper mo- 
tions) are given in Table 4. 

For our analysis we did not fit for the parallaxes of the 
background sources as they are primarily at distances of 
> 1 kpc (see Appendices A and B). Our final x 2 value 
for the fit was 228 for 50 degrees of freedom. The reduced 
X 2 deviates significantly from 1, indicating that we may 
be missing some source of error in our analysis. We have 
determined, though, that the majority of the excess x 2 
comes from four objects: the stars #21 and #201, and the 
galaxies #20 and #104. Without these sources, we obtain 
a x 2 value of 56 for 38 degrees of freedom, or a reduced 
X 2 of 1.5. This value is much more acceptable. The ques- 
tion, then, is why these four sources contributed so much 
to the x 2 ■ For the galaxies, they were not measured with 
the ePSF technique of AKOO, and are therefore subject 
to systematic errors not addressed here; we only include 
them to provide an "inertial" reference frame in the anal- 
ysis. The stars, #21 and #201, are among the brightest 
of the sources measured (Table 4) and may be among the 
closest (except for RX J1856.5— 3754, of course). As can 
be seen from their position residuals (Figure 2) they may 
have parallaxes of ~ 1.5 mas and therefore would not be 
expected to register well. They are also close to extremely 
bright stars (Figure 6) that may bias the position mea- 
surements. Together, these effects (non-zero parallax and 
mis-mcasuring positions) significantly increase the x 2 f° r 
the fit. Neither of these effects is important for the major- 
ity of the stars, which typically have residuals consistent 
with zero parallax (also see Appendix B). For the results of 
the registration, the inclusion or exclusion of these sources 
does not matter. 

We performed additional analyses to determine how ro- 
bust our measurements are, and included information from 
these analyses in the final estimates (see §2.4 and Ap- 
pendix A). 

The deviations of the scales and position angles from 
the nominal value were small but significant (see Table 1), 
unlike stated by W01. We find that the scale changed by 
w 0.03% from one epoch to another, and that the posi- 
tion angle changed by w 0?02. This is reasonable, given 
the fluctuations seen in other WFPC2 observations (due 
to thermal fluctuations in the detector and telescope as- 
sembly; AKOO). 

2.3. Absolute Astrometry 

Absolute astrometry was done relative to the USNO- 
A2.0 catalog (Monet 1998). We first determined centroids 
for all 571 USNO-A2.0 stars that overlapped with the av- 



The Parallax and Proper Motion of RX J1856. 5-3754 



5 



Table 3 

Distortion Corrected (x,y) Source Positions 



ID Epoch 





1996 Oct 




1999 Mar 




1999 Sep 






X 


y 


N' A 


X 


y 




X 


y 


A/a 




(pixels) 




(pixels) 




(pixels) 




i on 


269.63(3) 


173.83(2) 


o 


562.43(2) 


694.92(2) 


Q 
O 


304.53(2) 


170.38(2) 




102 


452.45(1) 


132.70(1) 


4 


379.55(3) 


732.33(1) 


8 


489.55(2) 


146.83(1) 


4 


103 


583.81(1) 


107.54(1) 


4 


247.36(2) 


755.08(2) 


8 


623.26(1) 


134.02(1) 


4 


104 


612.51(3) 


120.50(5) 


4 


218.50(5) 


741.63(7) 


8 


650.93(5) 


149.61(5) 


4 


105 


616.08(4) 


295.90(4) 


4 


218.82(4) 


565.99(4) 


8 


637.53(4) 


324.77(4) 


3 


106 


561.03(2) 


366.40(2) 


4 


275.145(9) 


496.71(1) 


8 


576.066(8) 


389.601(8) 


4 


107 


84.80(2) 


400.36(2) 


4 


751.80(2) 


471.926(9) 


8 


98.673(9) 


378.58(1) 


4 


108 


233.92(33) 


725.84(11) 


4 


609.04(8) 


144.08(9) 


7 


216.36(12) 


716.08(11) 


4 


110 


278.43(2) 


762.35(2) 


4 


565.18(2) 


106.41(2) 


8 


257.51(2) 


757.17(2) 


4 


111 


738.04(1) 


728.20(2) 


4 


104.91(1) 


131.54(1) 


8 


718.50(2) 


766.57(1) 


4 


112 


707.599(9) 


519.420(9) 


4 


131.642(9) 


340.811(9) 


8 


707.46(2) 


555.82(2) 


4 


113 


507.48(7) 


598.59(2) 


4 


332.92(2) 


265.71(3) 


8 


501.25(3) 


615.74(5) 


3 


114 


523.36(2) 


522.78(3) 


4 


315.58(1) 


341.14(1) 


8 


524.10(1) 


541.72(2) 


3 


116 


69.84(2) 


146.13(1) 


4 


761.95(2) 


726.29(1) 


8 


107.72(1) 


124.24(2) 


1 


117 


434.37(1) 


135.98(1) 


4 


397.66(3) 


729.44(2) 


8 


471.62(30) 


147.99(30) 


3 


118 


452.54(1) 


254.49(2) 


4 


381.66(2) 


610.43(2) 


8 


478.38(1) 


268.20(2) 


3 


119 


597.78(3) 


375.85(4) 


4 


238.33(3) 


486.51(3) 


8 


612.02(10) 


402.70(13) 


3 


127 


123.33(4) 


660.75(6) 


4 


718.28(4) 


210.87(4) 


8 


112.59(3) 


641.37(3) 


4 


128 


130.64(3) 


624.33(2) 


4 


710.32(2) 


247.15(2) 


7 


123.17(2) 


605.92(2) 


3 


129 


161.07(4) 


534.78(3) 


4 


677.93(7) 


335.92(3) 


8 


162.31(6) 


519.85(5) 


4 


201 


764.454(8) 


494.44(2) 


4 


74.173(8) 


364.652(8) 


8 


766.54(1) 


536.370(8) 


4 


J 


450.297(9) 


374.648(8) 


4 


385.946(8) 


490.481(8) 


8 


465.10(1) 


387.523(8) 


4 


19 


88.005(8) 


303.943(9) 


4 


746.47(2) 


568.192(8) 


8 


111.311(9) 


282.994(8) 


4 


20 


243.62(4) 


455.03(3) 


4 


593.76(3) 


414.32(3) 


8 


252.08(4) 


447.98(2) 


4 


21 


235.02(1) 


380.065(9) 


4 


601.15(1) 


488.972(9) 


8 


250.279(9) 


372.88(1) 


4 


23 


317.345(9) 


444.296(9) 


4 


520.60(2) 


423.45(1) 


8 


325.671(9) 


444.27(3) 


4 


25 


602.610(8) 


650.40(2) 


4 


239.092(8) 


211.774(8) 


8 


590.571(8) 


676.538(8) 


4 


26 


698.28(2) 


682.24(1) 


4 


144.40(1) 


177.714(8) 


8 


682.42(1) 


717.680(8) 


4 


28 


375.731(8) 


290.699(8) 


4 


458.94(1) 


575.797(9) 


8 


398.67(2) 


297.003(8) 


4 


X 


357.93(4) 


516.84(4) 


4 


472.30(4) 


333.92(4) 


8 


368.42(4) 


540.08(4) 


4 



a The number of independent measurements used to determine the mean position; see §2.1. 
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Table 4 

Reference Sources for Astrometry 



ID a 


a b -18 h 56 m 


<5 b +37°54' 


A/fqo6W 


Aa c 


AS C 


Ma 


^6 




(s) 


(arcsec) 


(mag) 


(arcsec) 


(mas yr ) 


100 


34.28 


-29.2 


25.0 


-12.9651(9) 


2.4976(9) 


0.1(5) 


5.6(6) 


102 


34.60 


-21.6 


24.5 


-9.1416(8) 


10.0946(9) 


-5.6(3) 


-6.6(3) 


103 


34.85 


-16.2 


23.9 


-6.2193(4) 


15.4760(3) 


-2.7(2) 


1.9(2) 


104 d 


34.96 


-15.6 


25.2 


-4.924(1) 


16.117(1) 






105 


35.49 


-20.6 


26.0 


1.3560(2) 


11.1550(2) 


-1.3(7) 


-0.9(7) 


106 


35.56 


-24.5 


22.4 


2.2431(3) 


7.1809(3) 


-2.7(3) 


0.9(3) 


107 


34.50 


-42.3 


23.7 


-10.3348(4) 


-10.5884(4) 


-1.8(3) 


-1.2(3) 


108 e 


35.83 


-46.5 


25.0 










110 


36.05 


-45.9 


25.1 


o.Uzoo^ / ) 


1 A r > r >QC\f'7\ 

— 14.zzoU(^ ( ) 




9 A ( A\ 
Z.4^4 J 


111 


37.07 


-28.7 


24.6 


20.1264(4) 


2.9624(5) 


1.5(3) 


4.8(3) 


112 


36.38 


-23.8 


23.6 


11.8772(5) 


7.9134(5) 


-2.4(2) 


-1.0(2) 


113 


36.12 


-33.1 


23.8 


8.8837(3) 


-1.411(1) 


0.1(9) 


6(1) 


11/1 
114 


Q/1 


A 


9/1 
Z4.Z 


6.6715(4) 


1.3405(4) 


-0.3(5) 


4.8(4) 


116 


33.71 


-35.5 


24.0 


-19.7242(6) 


-3.7665(6) 


-0.9(3) 


-3.3(3) 


117 c 


34.57 


-22.4 


25.8 










118 


34.97 


-25.1 


24.6 


-4.8365(6) 


6.5790(6) 


-1.9(3) 


-4.9(3) 


119 c 


35.68 


-23.6 


25.7 










127 


35.37 


-48.5 


26.3 


-0.043(1) 


-16.758(1) 


-1(1) 


1.1(9) 


128 


35.28 


-47.2 


24.9 


-1.1148(9) 


-15.4528(9) 


-1.2(5) 


-1.2(5) 


129 


35.09 


-43.5 


26.3 


-3.374(3) 


-11.782(2) 


4.5(8) 


1.8(9) 


201 


36.44 


-21.1 


20.4 


12.6471(3) 


10.6446(3) 


-0.7(2) 


0.6(2) 


J 


35.32 


-28.7 


21.0 


-0.6627(4) 


3.0372(4) 


-0.4(2) 


-0.3(2) 


19 


34.22 


-39.4 


21.9 


-13.6277(4) 


-7.6757(4) 


2.4(2) 


2.4(2) 


20 d 


35.05 


-38.3 


23.7 


-3.8055(8) 


-6.5591(8) 






21 


34.81 


-36.4 


21.0 


-6.6871(4) 


-4.7115(4) 


3.9(2) 


-2.3(2) 


23 


35.20 


-35.4 


23.8 


-2.0685(5) 


-3.6799(5) 


-5.0(3) 


-6.1(2) 


25 


36.51 


-31.3 


22.0 


13.4631(2) 


0.4201(2) 


0.2(3) 


-2.6(3) 


26 


36.84 


-28.9 


21.5 


17.3585(4) 


2.8446(5) 


2.2(2) 


-11.1(3) 


28 


34.89 


-28.9 


21.7 


-5.7783(5) 


2.8325(5) 


-0.4(2) 


-2.3(2) 


X f 


35.60 


-36.2 


26.1 


2.643(2) 


-4.504(2) 







a ID's are as in W01 where possible; source 201 was not present in W01; X is RX J1856.5— 3754. 
b Positions are measured at equinox J2000, epoch 1999.26. 

c Position offsets at epoch 1999.7, relative to the pointing center, for which our absolute astrom- 
etry yields a = -18 h 56 m 35=374 and 8 = -37°54'3l'.'71; see §2.3. 

d Extended source, probably a galaxy. 

e Rejected from the analysis due to poor position measurements; see Table 3. 
f For proper motions, see Table 5. 
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erage R-band image obtained in 2000 using FORS2 on 
UT#2 (Kueyen) at the Very Large Telescope (see van 
Kcrkwijk & Kulkarni 2001b; the image is composed of 29 
exposures of 135 s). We rejected 63 objects that were badly 
overexposed or had widths inconsistent with them being 
stellar. Next, the pixel coordinates were corrected for in- 
strumental distortion using a cubic radial distortion func- 
tion provided to us by T. Szeifert and W. Scifcrt (1999, 
private communication). Finally, the zero point position, 
the plate scale, and the position angle on the sky were 
determined, rejecting iteratively a further 87 objects for 
which the residuals to the solution were larger than 0'.'6 
(inspection of the images showed that virtually all of these 
were visual doubles, which are blended on the sky survey 
plates on which the USNO-A2.0 coordinates are based). 
For the 421 stars that pass our cuts, the inferred single- 
star measurement errors are 0'.'18 in each coordinate, which 
is line with the uncertainties expected for the USNO-A2.0 
catalogue (Deutsch 1999). Thus, we conclude that our as- 
trometry is tied to the USNO-A2.0 system at about O'.'Ol 
accuracy. 

We used the solution to determine the positions in the 
VLT R-band image of 19 stars from Table 4 (for the re- 
maining 7 - objects 105, 106, 112, 113, 118, 127, and 128 
- it was not possible to determine accurate positions, ei- 
ther because they were too faint or because they were too 
close to brighter stars). Using these epoch 2000.3 posi- 
tions and the fiducial epoch 1999.7 positions derived from 
the registration of the HST images above, we derived the 
pointing center for our reference HST image (we corrected 
for the difference in epoch using the observed proper mo- 
tions). We solved for zero point offsets only, i.e., the scale 
and orientation were held fixed to the values listed in the 
header. The inferred single-star measurement errors are 
0'.'025 in each coordinate, and the zero-point should thus 
be tied to the R-band image to better than O'.'Ol. 4 With 
this pointing center, we determined the absolute positions 
for all stars listed in Table 4. These should be on the 
USNO-A2.0 system to about 0702, and on the Interna- 
tional Celestial Reference System to about 0'.'2. 

2.4. Determination of Parallax and Proper Motion of 
RX J 1856.5- 3754 

With the three epochs registered, we compared the po- 
sitions of RX J1856. 5-3754 in each (see Table 3). We 
combined the initial estimates of the position uncertain- 
ties in quadrature with the uncertainties introduced by the 
registration. We fit for the proper motion and parallax 
of RX J1856.5— 3754 using a linear least-squares solution. 
The locations along the parallactic ellipse at each epoch 
were determined using the JPL DE200 ephemeris. 

As noted in §2.1 and §2.2, our limited number of mea- 
surements means that the individual position uncertainties 
have limited accuracy. Because of this, the uncertainty on 
the parallax derived from strict statistical considerations 
(1.7 mas) may not be correct. We have therefore estimated 
the parallax uncertainty using a variety of techniques (see 
Appendix A); these techniques have an rms of 0.4 mas, 
which we add in quadrature to find an rms of 1.8 mas. 

4 Leaving the position angle and scale free, we find changes of 0?09 
and 0.07%, respectively. The corresponding change in inferred point- 
ing center is ~ 07002. 



To be conservative we round this up, finding the overall 
l-o" uncertainty to be 2 mas, similar to the value found by 
W01. 

We present the results of the fitting for n, n ai and ^$ 
in Table 5. We also present the values for the derived pa- 
rameters of distance D and transverse velocity V±. The 
best-fit parallax is 7± 2 mas. We can exclude a null-result 
for the parallax at the 10~ 4 level. Our results are inconsis- 
tent with those of W01 at the 99.8% level. However, our 
best-fit values for the proper motions are entirely consis- 
tent with those presented in W01, and also agree with the 
orientation of the Ha nebula (van Kerkwijk & Kulkarni 
2001a). 

The severe inconsistency between the values of the par- 
allax derived by W01 and that derived by us is puzzling. 
The most obvious explanation for the discrepancy between 
our measured parallax and that of W01 is the significant 
difference in the way the data were processed: we used 
an ePSF that accounts for pixel phase errors and fit the 
data without manipulations such as rebinning or resam- 
pling. We also use more accurate distortion corrections, 
account for small changes in the scale and position angle 
of the observations from their nominal values, and account 
for the proper motion of the background objects. W01, on 
the other hand, first resampled the data, then shifted, then 
rebinned the data, and measured the positions with a PSF 
that is independent of pixel phase. Such analysis is liable 
to introduce even more pixel phase errors than those that 
were originally present. However, even if we follow the gen- 
eral method of W01 (resampling and rebinning, Gaussian 
fitting, old distortion corrections) we cannot reproduce a 
parallax of 16.5 mas (see Appendix A). 

The surest way to resolve the differences between our 
analysis and that of W01 will be with the 2001 March 
HSTd&ta, (not yet publicly available). This should allow 
further refinement of the proper motion and a more robust 
measurement of the parallax, with a final uncertainty of 
ss 1.5-2 mas. 



337 - 




330 - 



2 3 4 5 6 7 8 9 10 11 12 13 
71 (mas) 

Fig. 3. — Joint confidence contours for the parallax n 
(mas) and the magnitude of the proper motion fi (mas yr -1 ) for 
RX J1856.5— 3754, illustrating the covariance between these param- 
eters. Plotted are 68%, 90%, 95%, and 99% confidence contours, 
with the best-fit values indicated by the star. Note that the con- 
tours are for the joint confidence — the 1-D 68%confidence levels 
are given in Table 5. Compare to Fig. 2 from W01. 
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Table 5 
Motion of RX J1856. 5-3754 



Parameter 


Best-fit Values 


a a 


2.6435 ±0.0042 


5 a 


-4.5050 ±0.0030 


Ha (mas yr _1 ) 


328 ± 1 


M<5 (mas yr" 1 ) 


-58 ± 1 


7r (mas) 


7±2 


D (pc) 


140 ± 40 


H (mas yr 1 ) 


333 ± 1 


PA (deg) 


100.2 ±0.2 


V± (km s" 1 ) 


220 ±60 



a Fiducial positions at epoch 1999.7, 
relative to (x , y ) offsets given in Ta- 
ble 1. 



Note. — Best-fit values for /j, a , n&, 
it are determined directly from a lin- 
ear least-squares solution, and errors 
are l-cr/68% confidence. Best-fit val- 
ues and errors for the other parame- 
ters are derived from those for fj, a , [/,$ 
and 7r. 



As noted by W01, due to the small angle between the 
proper motion and the major axis of the parallactic ellipse, 
there is significant anti-correlation between the parallax 
and the magnitude of the proper motion. This is shown in 
Figure 3. We stress, though, that even with the latitude 
given by this anti-correlation we cannot accommodate a 
parallax of 16.5 mas. 

The position offsets with the proper motion subtracted 
are shown in Figure 4. The offsets are consistent with the 
best-fit parallax. The correlation between motion due to 
parallax and due to proper motion is also illustrated in 
Figure 4 — the proper motion direction differs from the 
position angle of the parallactic ellipse by only 20°. 

3. DISCUSSION 
3.1. Mass & Radius 

The most immediate impact of a revised distance for 
RX J1856.5— 3754 is in the interpretation of its spectrum. 
Pons ct al. (2002) used spectral fits to X-ray and broad- 
band data to determine a best-fit mass and radius, taking 
the previously published distance of d = 61 pc (W01) to 
convert the measured angular size Roo/d to a radius Roo, 
where Roo is the "radiation" radius determined from black- 
body fitting. The best-fit radius for all atmosphere choices 
was Roo = 6 km; causality then requires that the mass 
be less than this (in geometric units), giving M < 1M Q 
(Hacnsel 2001). These values are inconsistent with all neu- 
tron star equations-of-state (EOS), most of which have 
radii R > 10 km (Thorsett & Chakrabarty 1999; Lattimer 
& Prakash 2000). 

Our new distance pushes RX J1856.5— 3754 further 



6 
4 

2 

Sf 

CO 

E. 

< -2 
-4 




15 



10 



5 

A a (mas) 



Fig. 4. — Parallactic ellipse for RX J1856. 5-3754, showing the 
measured positions (points with error-bars) and expected positions 
(filled circles) at each epoch (after subtracting the best-fit proper 
motion). The inset shows an arrow that indicates the direction of 
the proper motion. The solid ellipse is for the best-fit parallax of 
7 mas, while the dotted ellipses are for the ±1-<T values of 5 and 
9 mas. This figure can be compared to Fig. 4 from W01. 



away, and therefore allows for larger radii and masses. For 
example, the best-fit uniform temperature blackbody fit 
(a uniform temperature is preferred by the limits on X- 
ray pulsations; Burwitz et al. 2001Ransom, Gaensler, & 
Slane 2001) has R^/d = 0.11 ±0.01 km pc -1 (Pons ct al. 
2002), which had implied Roo w 6.7 km for a distance of 
61 pc. Using our revised parallax value changes the radius 
to Roo w 15 km, giving a physical radius of R xs 12 km 
(for the canonical neutron star mass of 1.4M ; Thorsett 
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& Chakrabarty 1999). This is much more in line with the 
likely values for R determined by equations of state (e.g., 
Lattimer & Prakash 2000). 

Paczyhski (2001) has predicted that the passage of 
RX J1856.5-3754 near star #115 (from W01) in 2003 
will cause the apparent position of star #115 to change 
by w 0.6 mas due to gravitational lensing. Paczynski 
(2001) goes on to suggest that precise measurement of 
this deflection, perhaps by the new Advanced Camera 
for Surveys (ACS) aboard HST, could be used to mea- 
sure the mass of RX J1856.5— 3754 and provide an in- 
dependent estimate to constrain the equation of state. 
While we do note include star #115 in our analysis be- 
cause it is fainter than our detection threshold, we were 
able to estimate its proper motion. As expected from its 
faintncss, the measurement is not very accurate; we find 
H a = 1.2 ± 1.6 mas yr _1 , = 3.4 ± 1.6 mas yr _1 . The 
closest approach of RX J1856.5-3754 to star #115 is 0'.'2 
and should occur around April 2004. With the revised 
distance to RX J1856.5— 3754, we find that the Einstein 
ring radius will be (^e ~ 9 mas (eq. 2 of Paczynski 2001), 
giving a displacement of w 0.4 mas. This is somewhat 
less than the displacement predicted by Paczynski (2001). 
Based on the current WFPC2 data, it appears impossible 
to measure (in a reasonable number of orbits) the position 
of the star with enough precision to perform the suggested 
measurements (even assuming a displacement of 0.6 mas) 
with the ACS. 

3.2. Origin & Age 

W01 used the observed proper motion of RX J1856.5— 3754 
to trace it back towards the Upper Sco OB associa- 
tion, a source of supernovae during the last few million 
years (de Geus, de Zeeuw, & Lub 1989). He found, for 
d = 61 pc, that RX J1856.5— 3754 came within 16 pc 
of the association 0.9 Myr ago. W01 notes that for an 
unconstrained distance, the radius of closest approach is 
minimized for a distance d — 130 pc, entirely consistent 
with our measurement. This then gives a closest-approach 
0.4 Myr ago (with a radial velocity of +30 km s _1 instead 
of -60 km s _1 ), making RX J1856. 5-3754 half the age 
given in W01. This revised age resolves the minor discrep- 
ancy in the cooling history of RX J1856.5— 3754 found by 
Yakovlev, Kaminker, & Gnedin (2001). 

3.3. Energetics & Nature 

Another area where the distance enters is in modeling 
of the Ha nebula that surrounds RX J1856.5— 3754 (van 
Kcrkwijk & Kulkarni 2001a). In most of the modeling, the 
distance enters linearly and the factor of <~ 2 difference 
that we find here will not significantly change the conclu- 
sions. However, there are a number of quantities that have 
steeper dependencies on the distance. We examine each of 
these. 

In their estimate of the minimum pulsar wind energy 
loss E in the bow-shock model, van Kerkwijk & Kulkarni 
(2001a) find that E oc d 3 . We are able therefore to revise 

the limit to E ~ 8 x 10 32 c?f 40 erg s _1 , where the distance 
to RX J1856. 5-3754 is 140di 40 pc. This impacts on the 
estimates of the inferred spin period P and magnetic field 



B, giving 5 P ~ 1.5 s and B <~ 1 x 10 13 G. The new radial 
velocity is reasonably consistent with the inclination angle 
of 60° ±15° determined for the Ha nebula. 

If RX J1856.5— 3754 were powered by accretion, van 
Kcrkwijk & Kulkarni (2001a) find an accretion rate M = 
3 x 10 9 dgQ 3 ' 5 g s _1 , where the distance d = 60d 60 pc. For 
^60 = 2.3, as we find here, this then implies an available 
accretion power of ~ 3 x 10 28 g s . This is now almost 
four orders of magnitude less than the observed bolometric 
luminosity of RX J1856.5— 3754, which is revised upwards 
to ~ 2 x 10 32 d 2 40 ergs -1 , further supporting the claim 
that accretion cannot power the source. 

The final model for the Ha nebula considered by van 
Kcrkwijk & Kulkarni (2001a) is an ionization nebula. 
Here, RX J1856.5— 3754 can still be a pulsar, but its energy 
loss E must be small enough that any bow-shock nebula 
is smaller than the observed nebula. This leads to the 
result E~ 2 x 10 34 d? 4 5 erg s 1 , a much less constraining 
value than that given in van Kerkwijk & Kulkarni (2001a). 
However, we note that the ionization model has become 
less likely. With a distance of 140 pc the shape of the ion- 
ization nebula cannot be easily reproduced (van Kerkwijk 
& Kulkarni 2001a). 

3.4. Local Density of Neutron Stars 

Once the emission characteristics and size of one neutron 
star are well determined, they can be used to calibrate 
other sources. As an example, we derive a relation between 
the optical magnitude, X-ray blackbody temperature, and 
distance for isolated neutron stars, and apply it to the 
two brightest isolated neutron stars and PSR B0656+14, 
a nearby radio pulsar. 

The optical emission from RX J1856.5— 3754 and an- 
other isolated neutron star, RX J0720.4— 3125 (Haberl 
et al. 1997; Kulkarni & van Kcrkwijk 1998), is very well de- 
scribed by the Rayleigh- Jeans tail of a blackbody, although 
at a level slightly above that inferred from blackbody fits 
to the X-ray data (van Kerkwijk & Kulkarni 2001b). In 
particular, /„ oc v 2 in the optical regime. If we assume 
that the surface compositions of isolated neutron stars are 
similar, then /„ oc R 2 kTv 2 /d 2 , where R is the neutron star 
radius, kT is the effective temperature of the surface, and 
d is the distance. If the neutron stars all have the same 
radii, we can write 

rfioo = 1.4A/j^lO (y - 25 - 7) mas/5 
V 57 ev 

= i. 4 ^J^l()( B - 25 - 3 > mag / 5 , (3) 

where d = lOOrfioo pc, and B and V are the observed 
optical magnitudes. Here we use the best-fit blackbody 
temperature; while this is not always the best-fit overall 
to the X-ray emission (Pons et al. 2002), it is a simple 
model and the dependence of g?i o on kT is rather weak, 
so it suffices. We list the implied distances for the three 
other neutron stars that have thermally-dominated X-ray 
emission and optical counterparts in Table 6. As a general 
result of our new parallax, we decrease the local density 
of isolated neutron stars by a factor of ~ 10 with respect 
to that inferred from W01. 

5 Eqns. 8 & 9 from van Kerkwijk & Kulkarni (2001a) contain a small 
error: both P and B should go as d~ 3 / 2 , not d 3 / 2 as written. 
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Table 6 

Implied Distances to Neutron Stars 



Name 


kT 


V 


^100 


Rcfs. 




(eV) 


(mag) 






RX J1856.5-3754 


57 


25.7 


1.4 


1,2, this work 


RX J0720.4-3125 


79 


26.6 a 


3.0 


3,4 


RX J1308.8+2127 


118 


28.3 


6.7 


5,6 


PSR B0656+14 b 


73 


27.3 


3.3 


7,8,9 



References. — 1 - Pons et al. 2002; 2 - van Kerkwijk 
k Kulkarni 2001b; 3 - Haberl et al. 1997; 4 - Kulkarni & 
van Kerkwijk 1998; 5 - Schwope et al. 1999; 6 - Kaplan et al. 
2002, in preparation; 7 - Pavlov, Welty, & Cordova 1997; 8 - 
Koptsevich et al. 2001; 9 - Zavlin et al. 2001, in preparation 

a B magnitude. 

b We have taken the temperature of the dominant black- 
body component, and extrapolated the V magnitude from the 
Rayleigh-Jeans tail found in the UV. 



Note. — Calibrated using the parallax of RX J1856. 5-3754; 
see §3, Eqn. 3. 



For the radio pulsar B0656+14 we get a distance of 
cZioo ~ 3.3, near the low end of the values estimated 
through other means (250-800 pc; Mignani, De Luca, & 
Caraveo 2000) but still plausible. We understand that 
W. Briskcn will soon have a VLBA measurement of the 
parallax, which should provide independent confirmation 
of our results (W. Brisken 2001, personal communication). 

Data are based on observations with the NASA/ESA 
Hubble Space Telescope, obtained from the data Archive 
at the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astron- 
omy, Inc., under NASA contract NAS 5-26555. D.L.K. 
holds a fellowship from the Fannie and John Hertz Foun- 
dation, and his research is supported by NSF and NASA. 
M.H.v.K. is supported by a fellowship from the Royal 
Netherlands Academy of Arts and Sciences. J. A. ac- 
knowledges support from HST grant GO-8153. We thank 
S. Kulkarni and D. Frail for valuable discussions. 

APPENDIX A 

DETAILS OF ANALYSIS 

To test the robustness of our analysis, we performed 
the same general analysis but with variations on the input 
data set and analysis method. These variations included 
combinations of: 

• Using a six-parameter linear transformation be- 
tween the epochs (instead of the standard four- 



parameter transformation involving a shift, scale, 
and rotation). 

• Rejecting the stars with the largest position uncer- 
tainties. 

• Rejecting the stars whose derived proper motions 
had the largest uncertainties. 

• Rejecting stars #21 and #201, and the galaxies 
#20 and #104. 

• Rejecting the stars more than 300 pixels (15") from 
RX J1856. 5-3754. 

• Rejecting up to 10 stars at random from the sample. 

• Using the F606W ePSF instead of the F555W ePSF 
(see §2.1) 

All of these analyses gave entirely consistent results with 
rms variance of 0.4 mas, showing that our parallax mea- 
surements are not biased by any particular data points. 
Comparison of these parallax determinations allows us to 
better estimate the uncertainty in the parallax. To the for- 
mal error determined from the least-squares fit (1.7 mas), 
we add (in quadrature) the 0.4 mas rms found above. 

As another test, we used the same algorithm to mea- 
sure the parallaxes of all of the other stellar sources in 
the HST images. As expected from their photometric dis- 
tances (Appendix B), there were very few convincing par- 
allax measurements for these sources (Figure 2). The mea- 
sured parallaxes had a mean of —0.3 mas and a variance 
of 1.2 mas. The variance in the measured parallax was 
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roughly independent of the brightness of the star, down 
to the brightness of RX J1856.5— 3754, and is reasonably 
consistent with our estimation of the uncertainty of the 
parallax of RX J1856.5— 3754. We therefore believe that 
a conservative estimate of the 1-cr parallax uncertainty to 
be 2 mas. 

Finally, we performed the same analysis but with the ini- 
tial astromctry done using more conventional Gaussian fit- 
ting and with older WFPC2 distortion coefficients (Holtz- 
man et al. 1995; Trauger et al. 1995), like W01. Again, the 
results were consistent with those found using the more ac- 
curate ePSF fitting. 

APPENDIX B 
BACKGROUND SOURCES 



a number of the stars in the field could be main-sequence 
or giant stars in the Sgr dwarf galaxy. 

In Figure 6, we show the proper motion vectors for the 
background sources determined from the fitting. Most 
of the motions are small, < 5 mas yr" 1 , with the ma- 
jority being w 2 mas yr _1 . This is consistent with the 
magnitude of the proper motion of the Sgr dwarf galaxy, 
250 km s _1 at a distance of 25 kpc. The two galaxies 
(which were forced to have zero proper motion in our anal- 
ysis) provide an absolute reference for these proper mo- 
tions. The net proper motion of the background sources is 

<~ 0.5 mas yr _1 . As this is less than the uncertainty in the 
measured proper motion for RX J1856.5— 3754, the motion 
of the background sources should not bias the parallax of 
RX J1856. 5-3754. 



Figure 5 shows the color-magnitude diagram for the 
background sources that have reliable VLT photometry 
(van Kerkwijk & Kulkarni 2001b, a) and the distances im- 
plied from main-sequence fitting. Almost all of the sources 
are bounded by main sequences with distances from 2- 
25 kpc, with a number at implied distances of 10-15 kpc. 
Alternately, a number of the sources are consistent with 
red-giant branch stars at a distance ~ 25 kpc. 



14r 
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Fig. 5. — Color-magnitude diagram for background sources. R 
vs. B — R, with a Ay = 0.4 mag reddening vector and a model 
main sequences at a distance of 2 kpc (dash-dotted line), a main 
sequence at a distance of 25 kpc (dotted line), and a giant branch 
at a distance of 25 kpc (dotted line, marked "LC=III") from Cox 
(2000, p. 388, p. 392). 
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The implied distances of some of these sources, up to 
25 kpc, places then at the edge of the Milky Way and 
s=s 2.5 kpc above the disk, making them unlikely to be 
part of the Milky Way. The sources may be, however, 
in the Sgr dwarf galaxy. The heliocentric distance of 
the Sgr dwarf galaxy is w 25 kpc, and it has a line-of- 

sight extent of ~ 8 kpc (Ibata et al. 1997). The region 
near RX J1856.5— 3754 is « 7.5° from the center of the 
Sgr dwarf galaxy (Ibata et al. 1997), plausibly within the 
solid angle subtended by the Sgr dwarf galaxy. Therefore 



Fig. 6. — Derived proper motions for the background stars and 
galaxies used to register the three epochs, ovcrlaycd on the 1999.7 
HST image. The lines indicate proper motions going from the stars 
to the circles. The sizes of the circles indicate the errors in the 
proper motions. The arrow in the upper right indicates proper mo- 
tion with a magnitude of 10 mas yr —1 ; next to it is a 5" scale bar. 
North is up, and East is to the left. The arrow from the position of 
RX J1856.5— 3754 indicates the direction but only 1/20 the magni- 
tude of the source's proper motion (no uncertainties are plotted for 
RX J1856.5— 3754 as its proper motion is scaled such that the uncer- 
tainties would be invisible). All sources are labeled with identifiers 
from Table 4. 
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Fig. 2. — Residual positions for background sources with best-fit proper motion removed; see Figure 4 for RX J1856.5— 3754. Position from 
1996 October has dot-dashed cross; position from 1999 March has dashed cross; position from 1999 September has solid cross. All sources 
are labeled with identifiers from Table 4 in the upper-left corners. The last plot in the bottom row shows a parallactic ellipse for n = 1 mas. 



